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FOREWORD 


As  part  of  the  on-going  research  program  in  "Decision  Control 
Models  in  Operations  Research,"  Mr.  Richard  Ehrhardt  has  investigated 
the  structure  of  an  inventory  model  with  stochastic  replenishment  de¬ 
livery  lead  times.  The  existing  literature  on  this  topic  has  been 
corrected  and  extended  to  encompass  myopic  optimal  policies  for  finite 
planning  horizons  and  optimal  stationary  (s,S)  policies  for  infinite 
planning  horizons.  Efficient  methods  for  computing  infinite-horizon 
(s,$)  policies  are  also  given.  An  adaptation  of  the  Power  Approximation 
(Technical  Report  H7)  is  found  to  provide  excellent  performance  relative 
to  optimal  policies.  Several  sections  of  this  report  parallel  similar 
findings  in  earlier  reports.  Other  related  reports  dealing  with  the 
research  program  are  listed  on  the  following  pages. 


Harvey  M.  Wagner 
Principal  Investigator 


Richard  Ehrhardt 
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-  Abstract  - 

A  stochastic  lead  time  inventory  model  is  analyzed  under  the 
assumptions  that  (1)  replenishment  orders  do  not  cross  in  time  and  (?) 
the  lead  time  distribution  for  a  given  order  is  independent  of  the  num¬ 
ber  and  sizes  of  outstanding  orders.  This  study  corrects  errors  in  the 
existing  literature  on  the  finite-horizon  version  of  the  model  and 
yields  an  intuitively  appealing  dynamic  program  that  is  nearly  iden¬ 
tical  to  one  that  would  apply  in  a  transformed  model  with  all  lead 
times  fixed  at  zero.  Hence,  many  results  that  have  been  derived  for 
fixed  lead  time  models  generalize  easily.  Conditions  for  the  opti¬ 
mality  of  (s,S)  policies  are  established  for  both  finite  and  infinite 
planning  horizons.  The  infinite-horizon  model  analysis  is  extended  by 
adapting  the  fixed  lead  time  results  for  the  efficient  computation  of 
optimal  and  approximately  optimal  '$,S)  policies. 
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(s,S)  POLICIES  FOR  A  DYNAMIC  INVENTORY  MODEL 
WITH  STOCHASTIC  LEAD  TIMES 


We  consider  a  periodic  review,  single-item  inventory  system  where 
unfilled  demand  is  backlogged,  demands  during  review  periods  are  inde¬ 
pendent,  and  the  lead  time  between  placement  and  delivery  of  an  order 
may  vary  randomly.  We  require  the  joint  distribution  of  lead  times  to 
have  the  properties  that  (1)  replenishment  orders  do  not  cross  in  time 
and  (?)  the  marginal  lead  time  distribution  of  each  order  is  indepen¬ 
dent  of  the  number  and  size  of  outstanding  orders.  These  assumptions 
could  be  appropriate  in  practice  when,  for  example,  only  a  single  sup¬ 
plier  is  used  and  when  the  stocking  organization  places  orders  that 
are  small  and  infrequent  from  the  supplier's  point  of  view.  Replenish¬ 
ment  costs  are  comprised  of  a  setup  cost  and  a  cost  that  is  linear  in 
the  amount  ordered;  holding  and  shortage  costs  are  incurred  in  each 
period  as  a  function  of  period-end  inventory.  Our  criterion  of  opti¬ 
mality  is  minimization  of  the  expected  discounted  cost  incurred  during 
a  planning  horizon  which  may  be  finite  or  infinite.  Alternatively, 
when  we  consider  an  undiscounted,  infinite-horizon  model,  our  criterion 
is  minimization  of  the  expected  cost  per  period. 

A  finite-horizon  model  of  this  system  was  analyzed  by  Kaplan  fa] 
under  the  additional  assumptions  of  stationarity  of  all  model  para¬ 
meters  and  continuously  distributed  demand.  T ho  principal  results  won 
(1)  that  optimal  policies  can  be  computed  using  a  dynamic  pmst  an  havin' 
only  a  scalar  state  variable,  representing  inventor/  on  hand  plus  on 
order  before  ordering  and  (?)  sufficient  conditions  can  tie  found  Utr 


2 


the  optimality  of  base  stock  policies  and  (s,S)  policies.  Although 
the  findings  in  [4]  represent  a  breakthrough  in  the  study  of  stochas¬ 
tic  lead  time  systems,  the  results  had  two  complicating  features  that 
are  not  present  in  fixed  lead  time  models.  First,  the  parameters  of 
the  dynamic  program  were  not  simply  related  to  the  marginal  lead  time 
distribution.  Second,  sufficient  conditions  were  not  found  for  the 
optimality  of  myopic  ordering  policies. 

In  this  paper  we  correct  two  technical  flaws  existing  in  [4], 
allowing  an  intuitively  appealing  analogy  with  a  zero-lead-time  model. 
Then  we  establish  conditions  for  the  optimality  of  myopic  base  stock 
policies  and  present  generalized  conditions  for  the  optimality  of 
(s,S)  policies.  We  also  extend  the  model  to  encompass  infinite  plan¬ 
ning  horizons  and  show  that  optimal  (s,S)  policies  exist  under  standard 
conditions  on  the  cost  functions.  Finally,  we  present  efficient  algo¬ 
rithms  for  computing  optimal  and  approximately  optimal  (s,S)  policies 
in  the  infinite-horizon  setting. 

1.  MODEL  SPECIFICATION 

Wo  initially  consider  a  finite  planning  horizon  of  N  periods, 

numbered  backwards  from  the  end  of  the  horizon;  that  is,  the  final 

period  is  given  number  1,  and  the  initial  period  is  given  number  N. 

Demands  in  successive  periods  are  independently,  but  not  necessarily 

identically,  distributed.  Specifically,  let  the  demand  in  period  n 

be  represented  by  the  random  variable  D  with  mean  p  and  cumula- 

n  n 

five  distribution  function  v  .  Also,  let  1- ■  • ,  j  •  i,  be  the 

^  1  *  J 

be  the  convolution  of  ,...,<K.  We  assume  complete  backlogging  of 
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unsatisfied  demand,  so  negative  inventory  levels  are  permitted.  Also, 
there  are  no  losses  from  the  system  other  than  through  demand  satisfactim 
Costs  in  different  periods  are  related  by  the  single  period  discount 
factor  a.  Let  cn(z)  be  the  cost  of  ordering  z  units  in  period  n, 
wi  th 

C„(Z)  -  KnH(z)  +  cnz  . 

where 

0,  z  -  0 

H(z)  = 

1 ,  z  >  0  . 

We  assume  that  both  the  setup  cost  and  the  linear  portion  of  the 
replenishment  cost  are  paid  upon  delivery  of  the  order.  This  assumption 
does  not  entail  a  loss  of  generality,  since  payment  of  either  portion 
at  the  time  of  ordering  can  be  described  via  scaling  K  or  c  by  the 
expected  value  of  c ^ ,  where  L  is  the  random  lead  time. 

Let  Ln(x)  represent  the  holding  and  shortage  costs  in  period  n, 
where  x  is  the  ending  inventory  level  in  that  period.  Also  wo  do<  inn 
the  function 

Ei^Jy  -  <Dn  y  ...  4  nti_. n 
^t„.i(y-u)d.;.n>n.j(u)  .  (1) 

We  specify  replenishment  lead  times  as  identically  distributed 
random  variables  which  can  take  on  values  from  zero  up  to  a  fixed  mu/i- 
"!-  Let  a  given  lead  time  be  represented  h y  the  random  variable  I  , 
having  the  probability  distribution 
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«.i  =  P(L=i},  i=0, . . .  ,m  . 


The  joint  distribution  of  lead  times  is  characterized  by  our  assumptions 
that  (1)  replenishment  orders  do  not  cross  in  time  and  (2)  the  lead  time 
of  an  order  is  independent  of  the  number  and  size  of  outstanding  orders. 

following  the  development  in  [4],  we  focus  on  a  single  ordering 
decision,  and  let  V  be  the  number  of  outstanding  orders  immediately 
after  the  current  order  and  before  deliveries  are  received.  Let 
U  (  V)  be  the  number  of  outstanding  orders  after  the  current  delivery. 
Now  if  V  is  less  than  U  in  a  given  period,  we  know  that  only  the 
oldest  orders  ( U- V  of  them)  must  have  been  delivered.  This  is  a 
consequence  of  our  assumption  that  orders  do  not  cross  in  time.  Fur¬ 
thermore,  our  second  lead  time  assumption  (the  lead  time  of  an  order 
is  independent  of  the  number  and  sizes  of  outstanding  orders)  implies 
that  the  dynamics  of  order  deliveries  are  specified  by  a  sequence  of 
m  +  1  non-negative  numbers  {pQ ,p^ , . . .  ,Pm} -  It  is  shown  in  [4]  that 


P{U=j|V=i}  = 


P y  J  =  0  or  j  <  i 


i-1 

1  -  E  Pk>  3  =  i  , 

k=0  K 


where 


E  Pk  =  1  • 

k=0  K 


The  analysis  in  [4]  is  conducted  entirely  in  terms  of  the  p  's, 

as  opposed  to  the  ?,,'s.  An  expression  relating  the  p.'s  to  the 

J  J 
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e.j's  is  given  on  p.  495  of  [4],  but  it  is  incorrect.  We  derive  the 
correct  relationship  by  noting  that  the  probability  2-  that  {i  =  i) 
is  given  by  the  product  of  probabilities  that  it  is  not  delivered  ir. 
each  of  i  del i very  epochs  (0,1 ,2,. . . ,i-l  periods  after  ordering) 
and  the  probability  that  a  delivery  is  made  i  periods  after  ordering. 
We  have 


Although  our  derivation  of  (2)  is  original,  we  note  that  it  also 
appears  in  Nahmias  [5,  p.  913].  Since  the  substantive  results  of  [4] 
do  not  depend  upon  the  relationship  between  the  p.  's  and  ■: ^ '  s , 
they  remain  valid.  We  will  show,  however,  that  the  interpretation  of 
the  dynamic  program  in  [4]  is  simplified  by  (2). 

The  final  aspect  of  model  specification  concerns  the  costs  which 
must  be  included  in  computing  optimal  policies.  We  include  all  costs 
that  are  incurred  during  periods  11  through  1,  plus  those  that  occur 
in  the  following  in  periods  due  to  orders  placed  during  the  planning 
horizon.  A  terminal  reward  (or  salvage  value)  is  also  applied  to  the 
inventory  level  at  the  end  of  the  horizon.  This  differs  from  [4], 
which  only  considers  costs  incurred  in  periods  N  through  1,  and 
sets  the  terminal  reward  arbitrarily  at  zero  for  all  terminal  states. 

We  regard  this  to  be  another  flaw  in  [4].  Also,  we  will  show  that  our 
change  in  cost  accounting  allows  us  to  derive  conditions  for  the 
optimality  of  myopic  base  stock  policies. 
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2.  FINITE  PLANNING  HORIZONS 


The  central  finding  in  [4]  is  that  all  policy-dependent  costs 
can  be  included  in  a  dynamic  program  that  has  inventory  on  hand  plus  on 
order  as  its  only  state  variable.  Let  h^Cx)  be  the  minimum  expected 
discounted  cost  when  x  is  the  inventory  on  hand  plus  on  order  imme¬ 
diately  before  ordering  in  period  n.  We  have 


m 


h  (x)  =  min  £  aV[K  H(y-x)  +  (y-x)c  +  g  (i,y)] 
n  y>x  li=0  n  n  n 

( 


+  aEVi(y-DnH  >(3) 


h0(x)  =  -c0x  £  a1?,. 


Recursion  (3)  differs  significantly  from  those  given  in  [4]  only  in  that 
our  cost  accounting  includes  additional  terms  at  the  end  of  the  horizon, 
as  described  above  at  the  end  of  Section  1.  Now  let 


m 

£ 

i=0 


aV 

1 


m 

kn(y)  =  £  a\gn(i,y)/3  , 

i=0  1  n 


and 


Fn(*)  :  hn(x)/6  . 


7 


Recursion  (3)  can  be  rewritten  in  terms  of  this  notation  as 


(x)  =  min\K  H(y-x)  +  (y-x)c.  +  k 
’  y-x  (  n 


(y)  +  aFF  , (y-D  ) 
n  n- 1  n  ( 


r0(x)  -  -c0x 


(4) 


Notice  that  (4)  is  of  the  same  form  as  a  recursion  for  a  zero-lead-time 
system  with  k  (•)  representing  the  single-period  expected  holding  and 
shortage  costs. 

Following  the  approach  of  Veinott  [8],  we  establish  sufficient 
conditions  for  the  optimality  of  myopic  base  stock  policies.  Consider 
the  zero-lead-time  analogy  for  recursion  (4),  and  let  F  (xjY)  be 
the  expected  discounted  cost  in  periods  n  through  0  when  following 
a  particular  ordering  policy  Y  and  xn  is  the  starting  inventory  in 
period  j.  Also  let  iy. ,  7=1,..., N}  be  the  sequence  of  inventory  on 
hand  after  ordering  and  before  demand.  Then  one  can  show  that 


1  (x  |  Y ) 
tr  n  1 


n 

X 

i  - 1 


•<n_1  F 


K .  H(y . -x . )  +  G  •  (y . ) 


rn-1 

r 

i=0 


n-i 

<<  c,p 


i  i  + 1 


c  x 

n  n 


(5) 


where 


Gn(y)  =  (cR  -  <tcn  l  )y  +  kn(y)  .  (6) 

The  functions  Gn(*)  are  composites  of  expected  holding  and  shortage 
costs  and  the  linear  purchase  costs.  They  can  be  interpreted  (Veinott 
and  Wagner  [If)])  as  the  conditional  expected  holding  and  shortage  cost 
functions  of  an  equivalent  model  with  unit  purchase  costs  (  set  equal 
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to  zero.  We  shall  use  this  interpretation,  and  hereafter  refer  to 

f,  ( • )  as  a  conditional  expected  holding  and  shortage  cost  function. 

Consider  the  case  of  K  =  0  for  all  n.  It  follows  that  it 

n 

-G  (y)  is  uni modal  with  a  minimum  at  y  ,  and  if 
n  Jn 

y  -  D  -r  y  , ,  n=2 , . . .  ,N 

Jn  n  Jn-1 

with  probability  one,  then  it  is  optimal  to  order  max(yn  -  x  ,  0)  in 
period  n.  For  alternative  conditions  that  ensure  the  optimality  of 
myopic  base  stock  policies,  see  [8], 

An  interesting  parallel  to  the  fixed  lead  time  problem  arises  in 

the  solution  for  the  base-stock  values  {y  •  ,  i  =  l . N}  when  demand  has 

a  density  and  the  single-period  holding  and  shortage  costs  are  givpn  by 


Ln(x)  =  h  max(x,0)  +  p  max(-x,0),  n--m+l,...,N 


(7) 


Then  one  can  show  that  yn  is  a  solution  to 


<i„(y„)  ■  (p-Vacn_,)/(h+p)  , 


wtiere 


m 

o  (y)  =  V  a1?. .<l>  •  (y )/g  . 

n  pQ  i  n,n-i 

Notice  that  the  functions  a  (y)  are  linear  combinations  of  convolutt'd 

n 

demand  distributions  and  are  legitimate  distribution  functions  in  them 
own  right. 

When  the  model  does  not  possess  an  optimal  myopic  base-stock 
policy,  we  consider  the  function 


9 

=  Fn(xnlV)  ‘  [g  """’ViO  '  Vn_ 

=  £  <icn_1E  [KiH(yf-xf )  +  G.(y.)  . 

Notice  that  all  policy-dependent  costs  are  included  in  f^xjY). 
Therefore,  an  optimal  policy  can  be  found  by  computing 

fn(x)  =  min  f n ( x | Y ) 

using  The  dynamic  programming  recursion 

f  (x)  -  min^K  H(y-x)  +  Gn(y)  +  txEfn_1  (y-H  ) }■ »  n=l,...,N 
y  x 

•  (h 

f0(x)  =  o  . 

Lxpression  (8)  is  easily  recognized  as  a  standard  form  in  inventory 

theory.  Therefore,  conclusions  about  the  structure  of  optimal  polities 

are  immediate.  For  example,  if  =  0  for  all  n,  one  can  show  [t'J 

that  a  base-stock  policy  is  optimal  when  G  (y)  is  convex  for  all  n. 

★ 

the  base  stock  levels  are  given  by  the  values  ,  n~l . N}  that 

minimize  the  expression  in  braces  on  the  right-hand  side  of  . 

For  models  having  f<n  >  0  for  at  least  one  value  of  n,  we 
cite  two  theorems  from  Denardo  [1]  which  guarantee  the  optimality  of 
(s,‘>)  policies.  Let 

Jn(y)  -Gn(y)  *«Efn.,(y-<>„>  . 


: 

I 
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arid  define 


and 


Theorem  1 . 


Theorem  2. 


S  and  s  as  solutions  to 
n  n 


(y) 


sn  = 

inf{y  snlJn(y) 

=  K  +  J  (S  )} 
n  n'  n  ’ 

(>;!(•  th( 

■  following  three 

conditions  are 

(a) 

'!  (if) 

n  ‘ 

is  t -on  vex  when 

n  >  /,  and  -.7  j 

aont  inuous  and  K  -convex. 

(IT 

IS 

K  +  1 

■  aK  for  n='l , . 

(a) 

4  *  “ 

>  a  a  x  -*■  »  for 

n=ly  . . . ,N,  <m< 

t-f  (u- 

/'  J  is  finite ■ 
n+l  J 

for  n=l,...,N- 

77/.  'll 


K  +  j  ( U  )  y  X  j? 

n  n  n  >. 

Jn(x)  >  £>  a  , 


i.'ht  >‘f  ana'  r,  are  niner.  ha  (;‘)  tint!  ( I  nK 

n  n 

The  i  uitionr,  of  Theorem  1  remain  eai !  i  if  its  at 

iit  r‘ej> !•■:  ■<••  I  the  followin'.!  three  ■  •una’i  tionr: 

(a)  There  exist  a  '  a.,  ...  ■  rnoh  that  (_ 

is  Hon- irercusi.na  for  y  a  ana'  non-  :<  areas  i r 

ion  •  ■  ;  .  "orrorer,  a  (  >)  is  unit  •»;;*»>«?,  n 

n  n  ' 

i  *<!(;)■  1  i r/ j  ;  f’/  ll. 

>  n  n  '  n  ■■  1 

n  xx1 
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,  ■  oKfi  for  v.=  l 

(■■>  (:—!>  ,  J  in  J'in'  to  ton  / . V-  / . 

*  f!  *  >1+1 

Other  sets  of  conditions  for  the  optimality  of  (s,S)  policies  exist. 
See,  for  example,  Veinott  [9]  and  Schal  [7], 

3.  INFINITE  PLANNING  HORIZONS 


We  consider  an  infinite  horizon  version  of  our  model  in  which  all 
data  are  stationary.  Our  notation  is  simplified  in  this  setting  by  sup¬ 
pressing  subscripts  that  denote  period  numbers  whenever  the  quantity  of 
interest  does  not  vary  with  time.  Hence,  recursion  (8)  becomes 


fn(x)  -  min-jKH(y-x)  +  G(y)  +  uFf^  (y-O)j ,  n  1 


fQ(x)  0 


where 


G(y)  =  ( 1  ->x)cy  +  k ( y )  , 


arid 


in_ 

F(y)  £  '<  Vq(i  ,y)/H  . 
i=0 


Recursion  (11)  is  just  like  one  for  a  fixed- .end- time  model. 


with  G(y)  roprosentinq  the  conditional  expected  holdinq  and  short. iqe 
costs.  Hence,  we  know  that  if  G(y)  is  convex,  a  stationary  (s,s) 
Policy  is  optimal  in  (11)  as  n  approaches  infinity.  The  conclusion 


is  supported  by  the  argument  in  Iglehnrt  [8j,  which  also  establish*” 


the  existence  of 
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f(x)  =  lim  fn(x)  . 

n  w*. 

Also  from  [3],  we  know  that  f(x)  satisfies  the  functional  equation 

f ( x )  =  min<jt<H(y-x)  +  G(y)  +  aEf(y-D)l  .  (13) 

The  only  difference  between  (13)  and  a  fixed-lead-time  model  is 
in  the  function  G(*).  In  fact,  G(*)  can  be  expressed  in  the  same 
form  as  the  conditional  expected  holding  and  shortage  cost  function  of 
a  f ixed- lead-time  model  with  a  transformed  demand  distribution,  lie  have 

G(y )  =  (l-o)cy  +  k(y) 

m  . 

=  ( 1  -ot)cy  +  Y  a  £.g(i  ,y)/3 
i=0  1 


m  •  ~co  ^  »  » 

=  (l-a)cy  +  Y  a1*..!  L(y-u)d4>  U  '(u)/P  , 

*1=0  1  JO 

*  i 

where  T  J  is  the  j-fold  convolution  of  the  demand  distribution 
Hence,  G(-)  can  be  expressed  in  the  form 


G(y)  =  (l-a)cy  + 


(14) 


where 

m  .  ir  /  ‘  *  1  \ 

'V  (x)  =  Y  oiV*  1  (*)/(<  •  (!'•) 

x  i=0  1 

notice  that  the  function  f  (•)  has  all  the  properties  of  a  distribution 
function.  Therefore,  we  call  f  (•)  the  discounted  lead  time  demand 

ot 


distribution. 
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When  demand  is  discrete  and  6(*)  is  convex,  expressions  (14)  .md 

(lb)  allow  a  simple  adaptation  of  the  efficient  Veinott-Wagner  procedure 

[10]  for  computing  optimal  (s,S)  policies.  In  [10]  a  fixed  lead  time  \ 

is  assumed,  and  the  conditional  holding  and  shortage  cost  function  is 

*  ( >  + 1 ) 

computed  using  the  convoluted  lead  time  demand  distribution  T  ' 

To  adapt  [10]  for  our  stochastic  lead  time  model,  we  merely  substitute 
'P  in  place  of  <I>  in  any  computation  related  to  G(*)-  Specifi¬ 

cally,  the  important  expressions  in  [9]  that  require  modification  are 
(20),  (21),  (22),  (23),  (26),  (27),  and  the  two  unnumbered  expressions 
immediately  preceding  (21). 

We  have  performed  computations  using  the  procedure  described 

above.  We  summarize  the  results  below,  in  Section  4. 

We  conclude  this  section  with  a  discussion  of  approximately 

optimal  (s,S)  policies  for  the  infinite  horizon  model.  We  have  just 

shown  how  to  compute  optimal  (s,S)  policies  by  modifying  a  fixed- 

lead-time  procedure.  Basically  the  same  kind  of  modification  can  be 

used  to  compute  approximately  optimal  (s,S)  policies  as  well. 

lor  example,  consider  the  common  assumptions  of  ■  <  =  1  . ; nd 

linear  holding  and  shortage  costs  [as  given  by  (7)].  This  model  was 

analyzed  by  Roberts  [6],  who  used  asymptotic  renewal  theory  to  charac- 

★  ★ 

ter i ze  the  limiting  behavior  of  an  optimal  policy  (s  ,S  )  as  the 
parameters  K  and  p  grow  large.  He  obtained  the  following  expres- 

★  ★  ★  *  A 

'•  ions  ter  optimal  policy  parameters  s  and  11  S  -  s  ,  as  P 
grows  large: 

D  =  v2Kp/h  +  o(D  )  , 

C  *.  *(\+1 )  *  * 

\  *  l  n- s  )d!  u(u)  =  D  /(1+p/h)  +  o(P  )  , 

Js 


whore  A  is  the  fixed  lead  time,  u  is  the  single  period  demand  mean. 

•Jf  ic  ★ 

and  o(0  )/D  converges  to  zero  as  D  becomes  infinite.  Those  ex¬ 
pressions  were  used  by  Ehrhardt  [2]  to  construct  an  approximately  opt  ima 
policy  (the  Power  Approximation)  that  is  easy  to  compute  and  requires  fu 
demand  information  only  the  mean  and  variance  of  demand.  Specifically, 

the  Power  Approximation  requires  the  single-period  demand  mean  and 

p  *(  ul ) 

variance,  u  and  o  ,  as  well  as  the  mean  and  variance  of  T 

}i  and  ^ . 

We  suggest  modifying  the  Power  Approximation  for  our  stochast  ic 

2 

lead  time  model  by  replacing  and  with  the  mean  and  variance 

2 

of  ’1^,  p  ,  and  o^.  Specifically,  expressions  (13)  -  (16)  in  [l1] 

2 

require  this  change.  In  computing  p£  and  of ,  we  note  that  ■<  ~  I 
implies  that  P  =  1  also.  Therefore 


\  •  $0',dTi<u> 

-  E«,rud/(fti>(U) 

i=0  J0 


Me  -  =  (EL  +  l)u  , 


\  =  iU?dVl(u)  • 


=  (EL+l)o2  +  u?Var(L) 
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O 

Notice  that  p  and  are  merely  the  mean  and  variance  of  demand 
durinq  (L+l)  periods.  The  use  of  (16)  and  (17)  in  place  of  u;  and 
is  a  familiar  heuristic  approach  for  modifying  a  fixed  lead  turn- 
policy.  Until  now,  however,  this  approach  has  not  been  theoretically 
justified  for  periodic  review  systems. 

We  assess  the  effectiveness  of  the  modified  Power  Approximation 
in  Section  4,  below,  where  it  is  compared  with  optimal  policies  for  a 
variety  of  parameter  settings. 

4.  NUMERICAL  RESULTS 

We  have  performed  computations  using  the  procedures  described 
above  for  infinite  horizon  problems.  In  this  section  we  consider  a 
set  of  twelve  inventory  items  under  a  variety  of  assumptions  about  the 
distribution  of  lead  time.  First,  we  show  how  optimal  expected  costs 
vary  with  the  variance  of  the  lead  time  distribution.  Then  we  compare 
the  performance  of  optimal  policies  with  that  of  the  Power  Approximation 
as  modified  by  (16)  and  (17). 

Consider  a  system  of  12  inventory  items,  each  having  a  negative 
binomial  demand  distribution  with  a  variance-to-mean  ratio  of  3. 

Mean  demand  p  has  three  values,  2,  4,  and  8.  Each  item  has 
linear  holding  and  shortage  costs  as  given  by  (7).  Since  the  total 
cost  function  is  linear  in  the  parameters  h,  p,  and  K,  the  valur- 
of  the  unit  holding  cost  is  a  redundant  parameter  which  is  set  at. 
unity.  The  unit  shortage  costs  are  4  and  9,  and  the  setup  cost 
values  ore  3?  and  64.  All  combinations  of  these  parameter  settings 
are  included,  yielding  12  items. 
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We  consider  the  four  lead  time  distributions  displayed  in 
Table  I.  Each  is  a  symmetrical  triangular  distribution  over  the  range 
[0,4],  with  a  mean  value  of  2.  The  variance  of  lead  time  ranges  from 
a  minimum  of  0  for  the  deterministic  case  to  a  maximum  of  2  for  the 
uniform  distribution.  We  also  list  the  coefficient  of  variation  y  of 
each  lead  time  distribution,  which  is  defined  as  the  ratio  of  the 
standard  deviation  to  the  mean. 


TABLE  I 

Lead  Time  Distributions 


Probability  Mass  K. . , 

i=0,. . . 

,4 

EL 

Var  L 

i 

Y 

i 

2 

3 

4 

0 

0 

1 

0 

0 

2 

0 

0 

i 

0 

.250 

.500 

.250 

0 

2 

.50 

.354  ! 

.0667 

.2333 

.4000 

.2333 

.0667 

2 

1.00 

.500  i 

i 

.2 

.2 

.2 

.2 

.2 

2 

2.00 

.707 

Table  II  shows  optimal  total  cost  per  period  as  a  function  of 
parameter  values  for  each  of  the  four  lead  time  distributions  given 
in  Table  I.  Notice  that  the  total  aggregate  cost  of  the  12  items, 
increases  monotonical ly  with  lead  time  variance.  The  largest  lead 
time  variance  yields  an  optimal  total  cost  of  327  for  the  12  items, 
17  higher  than  the  deterministic  lead  time  cost  of  280.  When  costs 
are  aggregated  by  parameter  values,  we  see  that  the  larger  lead  time 
variances  produce  slightly  larger  cost  increases  for  items  with  a 


1 
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penalty  cost  p  of  9  than  for  those  with  p  equal  to  4.  Slightly 
larger  cost  increases  are  also  displayed  for  items  with  a  setup  cost 
K  of  32  as  opposed  to  those  with  K  equal  to  64.  The  bulk  of  the 
cost  increase,  however,  can  be  attributed  to  items  with  the  largest 
value  of  mean  demand.  Notice  that  items  with  p  equal  to  8  shows  a 
24'  increase  in  total  cost  from  126  for  the  deterministic  lead  time 
system  to  156  for  the  high  lead  time  variance  system.  T he  corres¬ 
ponding  percentage  increase  for  items  with  p  equal  to  2  is  merely 
8  .  This  fact  is  not  surprising,  since  we  have  held  the  demand  varianr.e- 
to-mean  ratio  constant.  Therefore,  items  with  the  largest  mean  demand 
also  have  the  largest  variance  of  demand,  yielding  especially  large 
values  of  a  in  (17). 


TABLE  II 

Optimal  Total  Cost  per  Period  of  12  Items  for 
4  Different  Lead  Time  Distributions 


Lead  Time 

Mean  Variance 

Total 

Aggregate 

Cost 

Costs  Aggregated  by  Parameter  Value 

Penalty  Cost 

4  9 

Setup  Cost 

32  64 

Mean  Demand 

2  4  8 

64  90  126 

65  93  136 

66  96  143 

69  102  156 

2  0 

2  1/2 

2  1 

2  2 

280 

293 

i  306 

327 

129  150 

135  159 

140  166 

149  178 

124  156 

131  162 

137  168  j 

149  179  j 

149 


178 


149 


179  j  69  102  156  I 


Table  III  lists  percentage  increases  in  expected  total  cost  per 


period  when  Power  Approximation  policies  are  compared  with  optimal 
policies.  The  Power  Approximation  yields  costs  within  a  few  tenths  of 
a  percent  of  optimal  for  all  parameter  settings.  This  level  of  per¬ 
formance  is  comparable  with  the  data  in  [2],  where  only  deterministic 
lead  times  were  considered. 


TABLE  III 

Percentages  Above  Optimal  Total  Cost  per  Period  for 
12-Item  Systems  Under  Approximately  Optimal  Control 


Lead  Time 

Mean  Variance 

Total 

Aggregate 

Cost 

|  Costs  Aggregated  by  Parameter 

Value 

Penalty  Cost 

4  9 

Setup  Cost 

32  64 

Mean  Demand 

2  4  8 

2 

0 

0.1 

0.2 

0.1 

0.1 

1 

0.3 

0.2 

0.2 

0.1 

2 

1/2 

0.2 

0.2 

0.1 

0.1 

0.2 

0.3 

0.1 

0.2 

2 

1 

0.2 

0.3 

0.1 

0.2 

0.1 

0.2 

0.2 

0.1 

2 

2 

0.3 

0.3 

0.2 

0.4 

0.2 

i 

0.0 

0.3 

0.4 
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